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Abstract 


Computer simulations of plasmas are relevant nowadays, because it 
helps us understand physical processes taking place in the sun and 
other stellar objects. We developed a program called PCell which is 
intended for displaying the evolution of the magnetic field in a 2D con¬ 
vective plasma cell with perfect conducting walls for different station¬ 


ary plasma velocity fields. Applications of this program are presented. 

This software works interactively with the mouse and the users can 
create their own movies in MPEG format. The programs were written 
in Fortran and C. There are two versions of the program (GNUPLOT 
and OpenGL). GNUPLOT and OpenGL are used to display the simulation. 

1 Introduction 

The huge advance in computer technology makes it possible to simulate and 
visualize complex physical phenomena taking place in stellar objects, e.g., 
the convective plasma cells on the sun. 

To understand the dynamo mechanism, most researchers consider cosmic 
plasma with a stationary motion, which leads to an induction problem. The 
aim is to hnd stationary states as solutions of the induction equation. Elsasser 
[5] , Weiss |T3] and Parker [10] performed the first two-dimensional simulations 
using symmetric velocity fields. 
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Simulations by Weiss and Galloway include dynamical effects and are gen¬ 
eralized to three dimensions for the kinetic case (Weiss [13], Galloway et al. 
[7], Galloway and Weiss |8]). In recent years, dynamo models have been 
improved by considering fully dynamical solutions of the induction equation 
taking into account the coupled mass, momentum and energy relations for 
the plasma. Recently, Mininni et al. [5] have considered the Hall current 
into the dynamo model. 

We can simulate convection cells by choosing convective velocity helds, which 
in turn help us understand the behavior of granules, mesogranules in the 
photosphere, supergranules in the photosphere and chromosphere, and labo¬ 
ratory plasmas. We also can investigate the not-well-understood phenomena 
of reconnection in this way [T]. 

The PCell program [2| helps visualize the magnetic held evolution in different 
convective plasmas. This program runs on Linux or Unix. The program can 
be downloaded from 

http://cinespa.ucr.ac.cr/software/xpcell/ 

The capability to easily create movies in MPEG format is one of the main 
features of PCell. 


2 The Induction Equation 

Maxwell’s Equations determine electromagnetic helds behavior in a cosmic 
huid 


eV-E 

= P 

(1) 


dB 


V X E 

~~dt' 

(2) 

V-B 

= 0, 

(3) 


dB 


V X B 

= /iJ -P 

(4) 


where p is the charged particles density. The ehects produced by the temper¬ 
ature gradient and charged particles density huctuations had been neglected. 
The general expression for the current density J in an isotropic homogeneous 
medium is: 


J = (tE + (jv X B -P pv 


(5) 
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where the right-hand-side terms are the conduction, induction and convection 
currents respectively. The conductivity a is considered constant in the whole 
plasma and v is the velocity held that describes the plasma motion. 

We can simplify the above equations by comparing the orders of magnitude 
of the quantities involved. We represent orders of magnitude with square 
brackets. 

In cosmic plasmas the velocity of the charged particles (mechanical velocities) 
are much slower the electromagnetic held velocity (speed of light). Therefore, 
we have: 


-J = [/^] < 1) (6) 

which implies that the orders higher than (3 can be neglected (non relativistic 
plasma). 

From the Faraday equation (the second equation of[T]), we can obtain an 
quantity with dimensions of velocity: 

= KB], (7) 

which is the velocity associated to the electromagnetic processes and satishes 
the condition: 


< [-y]- (8) 

The last two equations combined give us an estimation of the rate of the 
electric energy to the magnetic energy: 


■ eE^ 


■ E^ ' 





< in 


(9) 


If equation (| 8 ]) is fullhl, the electric component pE of the magnetic force that 
exerts the magnetic field over the plasma is negligible compared with the 
magnetic component J x B. 

The displacement current in the Ampere-Maxwell equation (the last equation 
of[T]) is negligible when is compared to the conduction current. The rate of 
both is given by 7 : 




where ujei is the electromagnetic frequency. Taking L as a characteristic 
length of the electromagnetic and mechanical phenomena [5], we have 
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and from (fTOj) we obtain 


H = 


(11) 


M < ( 12 ) 

For the Earth’s core 7 < 10“^® and for stars 7 <C 1 [1]. 

The rate of the convection cnrrent to the condnction cnrrent has the same 
valne 7 . From the Ganss law for the electric held (hrst eqnation of[T]) we 
have that [p] = [eE/L], therefore 


■ pv ■ 


' ev ' 

.aE. 


.aL. 


W- 


(13) 


According to eqnations ffTOj) and flTSl) . eqnation ([5]) and the last eqnation of 
(II]) simplify to 


J = (tE + (TV X B, (14) 

and 

V X B = /iJ. (15) 

respectively. 

Combining the latter two eqnations we obtain 

V X B = /iJ = paE + paw x B. (16) 

To eliminate E, we apply the cnrl to eqnation ffT 6 H and make the snbstitntion 
of Faraday’s law ([T]), yielding 

9B 

/i(T-^ = paV X (v X B) — V X V X B. (17) 

Finally, if we change the second term of the right side with the help of the 
well known identity and nse the magnetic Ganss law (third eqnation of[T]) we 
arrive to the indnction eqnation 

5B 

— = V X (v X B) + r^V^B, 
where rj = 1/pa is the magnetic viscosity. 


(18) 
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3 Behavior of the Induction Equation 

The rate of the first term on the right side of the induction equation to the 
second one is given by Rm = Lv/rj^ where the adimensional quantity is called 
the magnetic Reynolds number, in analogy to the Reynolds number for non 
conducting fluids. The bigger the plasma characteristic lengths is, the bigger 
the magnetic Reynolds number. 

If the hrst term is much bigger than the second, equation ffTSjl can be written 
as 


— = r^V^B. (19) 

This is the diffusion equation, which describes the magnetic field decay in a 
diffusion characteristic time 




( 20 ) 


for a plasma with spherical symmetry. It is of the order of one second for 
a one centimeter radii copper sphere, 10^ years for the Earth’s nucleus and 
10^*^ years for the Sun. 

The order of the magnetic Reynolds number can be written as 



tOel 


In. 

To 


( 21 ) 


where Tq is the time associated to plasma mechanical motion. Therefore, for 
times much slower than the diffusion time the equation simplifies to 


_ = V X (v X B) (22) 

which states that the magnetic flux trough any closed curve that moves with 
the local velocity of the plasma remains constant in time, i.e., the magnetic 
held lines are dragged by the huid {frozen field lines). 

When Rm ^ 1 the transport of the held lines by the plasma dominates over 
the dihusion, but if Rm -C 1 the held decays very fast and the dynamo ehect 
cannot take place. The behavior generated from the interplay of both terms 
for magnetic Reynolds number values between these limits is very interesting 
and to explore it is the aim of this work. 
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4 The Vector Potential Function 


The induction equation (1T8|) can be simplified for the two dimensional case 
if it is written as a function of the vector potential. We take the magnetic 
field and the velocity field limited to the x — y plane, then the magnetic field 
is obtained from the one component vector potential A = Ak as follows 


B = V X A 



after the substitution, equation (lT8jl yields 



(23) 


— ^ -u ■ + rjV'^A. (24) 

ot 

We define the position and velocity variables as function of the characteristic 
parameters (maximum velocity U and maximum length L) as follows 


u = u'U, X = x'L (25) 

with the analogous definition for the y coordinate. The spatial derivatives 
are given by 


d Id 

dx L dx' 

(26) 

d2 1 Q2 

dx'^ dx'^ 

(27) 


After the substitution of these relations and defining a the characteristic time 
of the mechanical motion tq = L/U. For instance, if we take following veloc¬ 
ity values for granules, mesogranules and supergranules in the sun, 

V = 900, 60 and400 m/s and L = 1.4 x 10^, 7 x 10^ and 3 x 10^ km, respec¬ 
tively, yield tq = 26min, 1.35and0.87days (Foukal [6], Sturrock et al. [n]). 
This time tq measures the time it takes the plasma to go from the bottom 
to the top of the cell. With t = f'r, we obtain 

where the primes have been removed for clarity. 

This is the equation we solve under the kinematic condition, i.e., there is no 
reaction of the magnetic field on the plasma, leaving the velocity field time 
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independent. This approach is valid if the magnetic energy is small compared 
with the kinetic energy of the plasma, that is 


Stt/j, 


1 2 


(29) 


5 The Visualization Program 

5.1 Description of the Program 

Equation fl28|) is solved using a fourth order difference schema in a two di¬ 
mensional cell with perfect conducting upper and lower walls (the magnetic 
held lines remain always tide to them) and periodic conditions at the lateral 
walls, i. e., each cell is surrounded by similar cells. 

The velocity held is taken to be incompressible, which allow us to dehne a 
stream function from which it can be obtained, we chose the following stream 
function [13] (see Figure 2) 



with (0 < X < 1, 0 < p < 1), m is an adjustable parameter (0 < m < 1) 
that allow us to select diherent velocity helds. 

When m = 1 it describe a single eddy [T3|, as m decreases two new sym¬ 
metrical eddies emerge from each side compressing the original eddie, at one 
time there are three eddies and as m get closer to zero, the central eddie 
disappear and remain just two eddies rotating in opposite directions. 

5.2 The Visualization Program 

The program called pcell is compressed in the pcell.tar.gz hie. The 
program was originally written in Fortran. We translated them into C with 
the help of F2C. The XFORMS library, which is used to design the control panel 
is also required. The data generated by the program is processed by GNUPLOT 
(PCell) or OpenGL (XPCell) to produce the simulation. 
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When the program starts, it creates a window: the PCell Control Panel. The 
user can control all items on it by clicking. 

5.3 The Control Panel 

The parameters are adjusted interactively by the user with the mouse. The 
control panel has the following items (see Figure 1): 

• The velocity held parameter selector: The user chooses the velocity 
held selecting a value of the m parameter, (0 < m < 1). 

• The magnetic Reynolds number selector: The user selects the magnetic 
Reynolds number, (0 < Rm < 1000). The values of U and L are chosen 
equal to one. 

• The time window: Displays the running time in units of tq = L/U. 

• The maximum running time selector: The time at which the program 
stops running can be chosen with this button. 

This items can be selected in any order but they should be selected 
before any other item of the program. 

• The Draw button: This button starts the program that calculates the 
magnetic held at each time interval (these data is stored in the hies 
pcelll.dat, pcell2.dat and so on) and after a short time, when the 
generation of data is hnished, the display window opens (see Figure 3). 

• The display window: The evolution of the magnetic held lines are dis¬ 
played in this window. When a given velocity held is selected it is 
shown here too. 

• The Stream button: Displays the mechanical motion of plasma (see 
Figure 2). 

• The Movie button: The users can create their own movies and with 
XANIM or a MPEG player they can display the movies. The name of the 
created movie is convection.mpg. 

• The Help button: Gives the user a program guide. 

• The Exit button: To leave the session. If the user wants to remove all 
pcell*.dat then type rm -rf pcell*.dat on the prompt. 
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6 Applications 

There are some interesting applications that the user can explore. Among 
these are: 

• The evolution of the magnetic held lines as function of the magnetic 
Reynolds number. 

• The mechanism of magnetic held dissipation and the reconnection phe¬ 
nomena. 

• The evolution of the averaged magnetic density as function of the mag¬ 
netic Reynolds number. 

• The maximum averaged magnetic density as function of the magnetic 
Reynolds number. 

• The stationary state of the averaged magnetic density and the way it 
is reached as function of the magnetic Reynolds number. 

• The time it takes to reach the maximum and the stationary state of 
the averaged magnetic density as function of the magnetic Reynolds 
number. 

6.1 Evolution of the magnetic field 

The stream function and the magnetic held twist evolution can be visualized 
easily with PCell. In Figure 2, an illustration of this evolution can be seen. 
Reconnection of the magnetic held appears in Figure 2b. 

6.2 Analysis of the magnetic field evolution 

For high Reynolds numbers, the convection generates a magnetic held com¬ 
ponent parallel to the plasma motion. The held begins to amplify itself, 
becoming stronger at the lateral borders of the cell (see Fig. 3). As this 
happens the spatial scale, where the variation of the held occurs, decreases 
linearly producing an increase of the dihusion. When this resistive term be¬ 
comes of the same order as the convective term, the amplihcation of the held 
stops. Therefore, the averaged magnetic density reaches a maximum (see 
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Figure 4). The twisting of the held produces regions of high magnetic den¬ 
sity at the same time that the diffusive term becomes bigger and reconnection 
of the magnetic hied lines occur expelling the held from the central region 
of the eddy reaching the magnetic density a stationary value (see Figure 2f). 
This stationary value is independent of the velocity held, but the way it is 
reached depends on the velocity held (see Figure 4). 

For velocity helds dominated by one eddy, the magnetic energy in the central 
part of the cell is very small compared with the density at the borders. As 
the held is expelled from the central part by dihusion, the magnetic density 
does not change much so the stationary value comes from the held located 
at the borders. For velocity helds with a strong inhuence of lateral eddies 
an additional region of high magnetic density appears in the middle of the 
cell, where the central eddy is located (see Fig. 3a - d). The intensity of the 
held is of the same order as the one at the borders; the twisting of the held 
as it disappears on both sides of the central region of high density makes 
the average magnetic density oscillate until it reaches the stationary value. 
When the velocity helds are practically two eddies (small m), the magnetic 
held at the borders is not very strong and the generation of the central high 
density zone helps to reach the maximum value which is the same as the 
stationary one (see Figure 4). 

It is important that the average lifetime of the granules and mesogranules 
is very small compared to the characteristic time (~ 0.31 and O.OTtq, respec¬ 
tively). Therefore for these motions, the behavior shown would not evolve 
because the structures have disappeared long before. On the other hand, in 
supergranules (~ 23ro) the stationary state is reached. 

According to the above description for single eddies, a simple expression for 
the maximum averaged energy density and the Reynolds number can be ob¬ 
tained as where Bq is the initial averaged magnetic held 

(Weiss [T3], Carboni [T]). In Figure 5a, it can be seen that this behavior is 
still valid for plasmas with interacting eddies. The Reynolds number expo¬ 
nents obtained for the values of m = 1.0, 0.7, 0.3 and 0.1 are, respectively, 
0.59, 0.55, 0.54 and 0.67. In the same way, a relation between the station¬ 
ary and averaged energy density is obtained (Weiss [13], Carboni [T]) as 
B‘1^ ~ Figure 5b shows the result of the simulations which agree 

very well with this relation even for interacting eddies. The Reynolds num¬ 
ber exponents obtained in this case are 0.42, 0.42, 0.43 and 0.48 in the same 
order as above. 
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7 Conclusions and Future Work 

The kinematical dynamo model presented here prodnces intensification of 
the magnetic field by the indnction of plasma trying to move across field 
lines in regions of small eddies between big convective zones in addition to 
the usnal one accumulated on the regions between cells. Although this model 
does not include the feedback of the J x B force on the inducting motion 
and has no rotation profile, the generation of a toroidal field is clear. Similar 
oscillatory behavior of the magnetic field for bands of asymmetrical eddies 
has been found for groups of four-cell convections by Zegeling [15] working 
on the same principle described here. 

The model could be implemented to solve the full dynamo problem, which 
involves the simultaneous solution of the induction equation along with the 
equations of motion, continuity and thermodynamics. The inclusion of the 
J X B term produces freezing of the magnetic field to the electron flow instead 
of the bulk velocity field (Mininni et al. |H]). 

Implementations of the numerical methods will allow us to expand the study, 
including the effects of nonlinearity and chaotic motion at Reynolds numbers 
typical of astrophysical problems where self-organization emerges (Chang et 
al. [3], Valdivia et al. [I2])- New approaches have been developed to calcu¬ 
late, in a more efficient way, convective cells. One strategy is the adaptive 
grid method which, based on coordinate transformations between physical 
and computational coordinates, automatically tracks and spatially resolves 
nonlinear structures (Zegeling and Keppens |16jb 

Future Work. The program can be improved in the following way: 

• Including more velocity fields. 

• Adding the magnetic density averaged over the cell, represented as 
function of time. 

• Expanding to three dimensional cells and using other shapes like hexag¬ 
onal cells (this shape appears as stable patterns in some fluids). 

• Considering mechanical-electromagnetic interaction between the plas¬ 
ma and the field. 

• Exploring more complex behaviors such as Chaos. 

• Using the more advance program XANIM to produce a better animation. 
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Figure 1: PCell Control Panel 



































Figure 2: The stream function for m = 0.1 (a) and the correspond¬ 
ing time evolution of the magnetic field with Rm = 1000 (b-f) at times 
t = 1, 1.5, 2, 2.5, 3. 
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Figure 3: Magnetic field states for different magnetic Reynolds numbers and 
plasma velocity fields: (a) Rm = 200, m = 0.1, (b) Rm = 500, m = 0.3, (c) 
Rm = 800, m = 0.7 and (d) Rm = 1000, m = 0.1. 
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Figure 4: Square of magnetic field as a function of time for different m values 
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Figure 5: Stationary (a) and maximum magnetic energy density (b) as func¬ 
tion of the magnetic Reynolds number for different m values. 











